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We investigate the dynamical formation and evaporation of a spherically 

symmetric charged black hole. We study the self-consistent one loop order 
semiclassical back-reaction problem. To this end the mass-evaporation is mod- 
eled by an expectation value of the stress-energy tensor of a neutral massless 
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. scalar field, while the charge is not radiated away. We observe the formation 



of an initially non extremal black hole which tends toward the extremal black 
hole M = Q, emitting Hawking radiation. If also the discharge due to the 



. instability of vacuum to pair creation in strong electric fields occurs, then 

§; 

until the scale where the semiclassical approximation breaks down. We cal- 
culate the rates of the mass and the charge loss and estimate the life-time of 
q-i the decaying black holes. 



I. INTRODUCTION 

Extremal black holes (EBHs) play an important role in black hole thermodynamics. The 
vanishing of the surface gravity of such black holes implies zero Hawking-Bekenstein tem- 
perature. This suggests on the tight connection of EBH in black hole physics with zero 
temperature states in thermodynamics. The zero temperature EBH is, however, unattain- 
able in a finite sequence of physical processes QXJ] . This means that the dynamical evolution 
of a non extremal black hole towards the EBH state continues an infinite (advanced) time. 
The full non linear evolution towards an EBH presents a step towards understanding the sta- 
tus of EBHs in the black hole thermodynamics. The solution provides the entire dynamical 
spacetime which emerges in this case. 
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A generic black hole carries an angular momentum and an electric charge. The formation 
and evaporation of such a black hole is a challenging and extremely difficult problem both 
analytically and numerically. That is because the problem is necessary three dimensional. 
Therefore, we consider here a simpler spherically symmetric problem of a formation and 
a decay of a charged black hole. The EBHs in this simple model have all the features 
of a more general Kerr-Newman black hole. In earlier studies [@dll of the evaporation of 
Reissner-Nordstrom (RN) black holes an initially singular spacetime-the static extremal RN 
black hole-was considered. This extremal black hole was exited above the extremality by 
an infalling matter and then allowed to decay back toward the Q = M state by emitting 
Hawking radiation. Here we pose the problem of the evaporation of a RN black hole in 
somewhat different fashion. We formulate the problem of the formation and the evaporation 
of a spherically symmetric charged black hole beggining with the collapse of an initially 
regular self-gravitating charged matter distribution. We include quantum effects, which are 
considered here in the one loop order, in a self-consistent manner. To do so, we introduce the 
expectation value of a scalar field stress-energy tensor, which acts as an additional source 
in the semiclassical Einstein equations (see e.g. ||). We regard the charge as stable, not 
allowing it to be radiated away Q The semiclassical equations are integrated numerically 
to give the structure of the spacetime for this configuration. We find that the evaporation 
proceed to a stable end-point corresponding to the extremal, M = Q, charged black hole. 

In the next stage we relax the assumption of a stable charge, allowing the black hole to 
discharge via the Schwinger pair creation process in addition to Hawking evaporation. The 
characteristic quantity governing the Schwinger process is the critical pair-creating field, 
E c = Trm 2 c 3 /eh, with e and m are the charge and the mass of the produced particles. 
Suppose that the critical field is reached below the maximal radius of the outer apparent 



This can be realized in practice when the charged particles are sufficiently massive and there is 
no significant creation of them in the Hawking or Schwinger processes. 



2 



horizon (i.e. the apparent radius of a non-evaporating black hole), but before the EBH 
state is approached. In this arrangement the black hole evaporates, radiating away its 
mass, the outer apparent horizon shrinks until the electric field upon this horizon reaches 
the critical pair-creating value. Then electrically charged pairs are produced intensively. 
Particles having the same polarity as a black hole are repulsed to infinity. Particles with the 
opposite polarity are captured by the black hole reducing its charge. Therefore, instead of 
reaching the stable EBH state the black hole evaporates and discharges simultaneously. The 
black hole looses its mass and the charge, disappearing completely (until it reaches scales 
where the the semiclassical approximation breaks down). It is of interest to note that the 
life-time of such decaying black hole is longer than that of a neutral evaporating black hole 
of the same mass. 

The article is organized as follows. In section |I| we introduce the physical model of 
and write down the evolutionary equations. The characteristic problem and the numerical 
scheme are described in section [TTT]. Section [IV] is devoted to our numerical and analytical 
results and in section |VJ we summarize our findings. We use units in which c = G = h = 1. 

II. THE MODEL 

We introduce a simple spherically-symmetric four dimensional model which captures the 
essence of the realistic behavior of the system. The semiclassical back-reaction is included in 
one loop order by introducing the expectation value of the stress-energy tensor, which arises 
from the quantization of a 2D massless scalar field in 2D curved spacetime || . To mimic the 
AD radial dependence we divide this tensor by 4irr 2 . This stress-energy tensor offers a good 
approximation to the full four dimensional theory since (i) most of the energy is carried away 
in S— waves J7|, and (ii) the inclusion of higher angular modes does not change significantly 
the behavior derived from the S— wave approximation ||, see, however 0. Even though the 
tensor, which we utilize, arises in a 2D theory and not from the spherical reduction || of 
the full 4D theory, it gives a feeling of what the back-reaction would be, and, indeed, this 
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tensor engenders the evaporation of black holes (see e.g. Ref. RllOfl). 

Our final aim is the numerical integration of the evolution equations. To simplify the 
numerical procedure we apply several convenient tools. The charged collapsing matter is 
simulated by a complex massless scalar field. The "massless" nature of the fields makes it 
very convenient to use the double-null coordinates. The line element is taken to be of the 
form: 

ds 2 = —a{u,v) 2 dudv + r{u,v) 2 dVt 2 , (1) 

where dfl 2 is the unit two-sphere. There is a coordinates gauge freedom. For the time being 
u, v are just general light-cone coordinates which would be specified by fixing a gauge (see 



section |T|) 



We formulate the set of coupled Einstein-Maxwell- Klein-Gordon equations as a first order 
system of PDEs. The numerical integration of the first order system functions very well both 



for the neutral [II], and for the charged [p~2| , [T3[| situations. It is convenient to define the 



auxiliary variables: 



d=—, h = —, f = r u , g = r v , w = s u , z = s v , (2) 
a a 



wherein s is the complex scalar field divided by y/Aw. We have adopted the notation W x = 
dW/dx for partial derivatives of any function W = W(x,y). 

The Hawking radiation is modeled by the expectation value of the quantum stress-energy 
tensor of a 2D massless scalar field. In the light-cone coordinates (H) the only non vanishing 
components of this tensor (divided by Airr 2 ) are: 

P 

(T uu ) = ^(h u -h 2 ), (4) 
(T m ) = ^(d v - d 2 ) , (5) 



wherein P is a constant, which is proportional to the number of massless scalar fields. P 
controls the rate of the evaporation and it also defines the length scale where the semiclassical 
approximation breaks down. 



We write the coupled set of the semiclassical Einstein-Maxwell-KG equations: 
Einstein equations: 
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Maxwell equations: 



a 2 q 

M1 = a v - 2^f = ' 



M2 = q v 



ier 2 (s*z — sz*) = 



Where a(u, v) = A u is the electromagnetic potential and q{u,v) is the charge. The complex 
scalar field (Klein-Gordon) equations expand to: 



le 



SI = rz u + fz + gw + iearz + ieags + —a 2 qs = 

Ar 



S2 = rw v + gw + fz + iearz + ieags + —a 2 qs = 

Ar 



Finally, the definitions (0) are rewritten as: 
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These equations should be accompanied by a specification of the initial data and the suitable 
boundary conditions. 
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III. INITIAL SETUP AND THE NUMERICAL SCHEME 



The characteristic initial value problem for the complex scalar field is a straightforward 



generalization of that given in Ref. [13]. We have formulated it in [13] and describe it briefly 
here. 

We choose the initial characteristic surfaces to be: the ingoing v = const = Vi hyper- 
surface, and, the outgoing u = const = Ui hypersurface. The coordinate gauge freedom 
is fixed by constraining the r to be linear with u or v along the characteristic hypersur- 
faces. Namely, on u = Ui segment we choose g = r v = 1, on v = v$ segment we choose 
/ = r u = r u0 . To get r along initial surfaces it is necessary to supply r = r(ui,Vi) that 
serves as a free parameter. For convenience we choose Ui = 0, Vi = ro- Therefore, one 
obtains: r(wj, v) = v , r(u,Vi) = ur u0 + r . 

The domain of integration in the v direction confined in the compact region: vi < v < Vf. 
The origin r = 0, for simplicity, is not included in the integration domain. This is achieved 
by an appropriate choice of the final outgoing hypersurface u/. The scalar field distribution 
is chosen to be non- vanishing only in the compact segment v , < v < v p (v p = V2,v' 2 ) and is 
taken in a shape of a pulse, which is smoothly matched at the endpoints (vi and v p = v 2, v' 2 ) 
to the rest of the integration domain: 

. . r,/ V — Vl\ . . 2 / V — Vl \ , . 

s[Ui, v) = Asm 7r \+tBsm 7r— , (6) 

V v 2 — vlj \ v 2 — vlj 

where A, B are the amplitudes of the pulse. The initial values of integrated variables are: 

An . / v - v x \ iBn . ( v- v x 
— sin 2tt— H sin 2ir 



v' 2 — Vi \ V 2 — V 1 J V2 — V 1 V V2 — V \ 

w(u,Vi) = . (7) 

The initial value for d(ui, v) can be obtained by a numerical integration of the constraint 
equation E2 along the initial u = Ui hypersurface, together with the choice a(ui,Vi) = 1. 
We neglect the quantum stress-energy tensor on the initial hypersurfaces. This is justified 
as we verify later in the numerical solution. From E2 we obtains the initial value for d: 
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7/ \ A 2 7l 2 V nf V — Vi \ 



B 2 7T 2 V 



V — V i 



sin ^ 27T— . (8) 

With no sources the spacetime is Minkowski spacetime for v < t>j. The boundary value 
for a(u, is obtained from the second constraint equation £"3 to be: a(u, Vj) — 1. The rest 
of the boundary setup is q(u, Vi) = 0, a{u, Wj) = 0. 

In our coordinates the mass-function is: M(u, v) — | ( 1 + ^ + ^ r M r i))- Since it vanishes 
for Minkowski spacetime (in the region v < V{) one can calculate: r U Q = It should be 
noted that for v ^> M our ingoing null coordinate v is closely related (proportional) to the 
ingoing Eddington-Finkelstein null coordinate v e . The u coordinate measures the proper 
time of an observer at the origin \ T3 \. 

The numerical scheme used to evolve the initial data is similar to the one used for 
the classical equations [I3~|. At each step we evolve d and z using El and 51 from the 
hypersurface u to u + du. We calculate the derivative d v . Then we solve the appropriate 
equations for the rest of quantities along the outgoing null rays u + du = const, starting 
from the initial outgoing hypersurface v = Vi. We integrate equation Dl to find a, then 
we solve the coupled differential equations D2 and E2 to obtain r and g. Next, equations 
D3, M2, Ml are integrated to obtain s, q and a. Finally, the differential equations E3 and 
52 are solved for / and w, respectively. 



This integration scheme uses several different numerical methods [13] to evolve the initial 
data. To evolve the quantities in the w-direction we utilize the 5-th order Cash-Karp Runge- 
Kutta method. The differential equations in the v direction are solved using a 4-th order 
Runge-Kutta method. The integrations in the v direction are performed using a three-point 
Simpson method. The new feature of the current scheme relative to the classical one is the 
necessity of calculating of partial derivative d v in equation E2. This derivative is calculated 
using a Savitzky-Golay smoothing filter. 
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IV. RESULTS AND DISCUSSION 



A. The Evaporation 



As a test we check that our code reproduces the uncharged black hole evaporation ||,[T0 . 
To do this, we solve numerically the set of the coupled Klein-Gordon-Einstein equations. 
The complex scalar field is replaced by a real one. We set B = 0, e = 0, but A ^ 0. We set 
the free parameters to be: r = v i = 10, Vf = 200, Ui = 0,t>3 = 20, P = 0.01 and A = 0.65. 
In Figure [1| we show the dynamical spacetime obtained in this case. 





[a) (b) 
FIG. 1. The evaporation of a Schwarzschild black hole, (a) The radius as a function of v along 

a sequence of an outgoing u = const null rays. The rays are emitted from the bottom left corner 

at different u = const moments. In the late u times the rays begin to curve towards the r = 

singularity. There are, however, rays which curve back toward the singularity, but then managed 

to escape to infinity, indicating the decrease of mass - the evaporation, (b) The contour lines of 

the radius in uu-plane, in a vicinity of the event horizon (defined as a last ray to escape the r = 

singularity). The radius decreases from the bottom to the top. The apparent horizon occurs in 

the locus of r v = and its radius decreases - a clear evidence for evaporation. 



Figure |T](a) depicts the radius as a function of the outgoing null coordinate v along 
a sequence of an outgoing u = const null rays. The evaporation is mandatory from this 
picture. Rays which initially incline towards the central r = singularity manage at some 
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later moment (at the apparent horizon, where r v = 0) to escape to infinity. The evaporation 
is indicated by the decreasing radius of the apparent horizon. Figure |l](b) displays the 
radius contour lines in the w-plane. The apparent horizon is indicated by vanishing of the 
r v derivative. Again, the decrease of this radius is clear. In passing, we note that our null 
coordinates are not very convenient to resolve the evaporation of a neutral black holes. The 
evaporation is contained within a tiny u lapse. The convenient choice is the null coordinates 
defined at infinity ||. In these coordinates the evaporation takes an infinite (retarded) time. 
We did not used these coordinates since our target is the entire spacetime of an evaporating 
charged black hole. Therefore, we need a set of coordinates which are regular on the outer 
horizon in order to push the integration inside the black hole. Our original coordinates, 
being Kruskal-like, supply this set. We find out that our results are in agreement with the 
previously established ones ||[10| and turn to the main concern of this work - the evaporation 



of charged black holes. 

We have set the rest of the free parameters as: v% — 5, 1)% — 14, v' 2 — t>3, f / = 2000, A = 
0.35,5 = 0.43 and e = 0.151. The Figures below are obtained with this specific choice of 
the parameters. Figure ^|(a) displays the radius r(u, v) as a function of the advanced time 
v along a sequence of u = const outgoing null rays. The retarded time u increases from the 
rightmost ray to the leftmost one. The rays are emitted from the surface of a collapsing 
matter (at the left bottom corner). The first rays, which escapes to infinity, are unaffected 
by the matter indicating the asymptotic flatness of the spacetime. As the collapse proceeds 
the rays begin to curve back towards the origin. Without evaporation there would be a 
ray which does not escape to infinity, but stays at a constant radius, this ray indicates the 
event horizon. The rays emitted after that would asymptotically (v — > oo) become vertical 
indicating the formation of the Cauchy horizon (CH) |L2],[T3|. Finally, the rays emitted in 
late retarded time would inevitably fall into the strong central singularity. In Figure [5] we 
depict the rf-diagram obtained in the classical without evaporation. 

Inclusion of the evaporation changes this picture, specifically the inner and the outer 
horizons are not stationary null hypersurfaces. Figure |2](c) depicts the enlarged part of 



Figure ||(a) - the part where the outgoing rays curve to be swallowed by the formed black 
hole, but then managed to escape to infinity. The locus of the turn-points of these outgoing 
rays indicates the outer apparent horizon. Figure 0(b) displays another enlarged part of 
Figure 0(a)- the part where the outgoing rays tend towards the CH. The incline of these, 
classically vertical, rays towards the outer horizon (above the thick dashed line, see Figure) 
can be observed. The turn-points (i.e. vanishing of the r v derivative) of these rays indicate 
the inner apparent horizon. 
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FIG. 2. The evaporating charged black hole. Figures depict radius r(u,v) as a function of the 
advanced time v along a sequence of u = const outgoing null rays. The retarded time u grows from 
the rightmost ray to the leftmost one. (a) The entire spacetime for a charged collapse, (c) The 
part of (a) in the vicinity of the outer apparent horizon (signaled by vanishing of the r v derivative), 
which displayed as a thick dashed line. The horizon initially grows, indicating absorption of a 
matter and formation of a black hole. After that the horizon contracts, indicating the loss of mass 
by the evaporation, (b) The part of (a) in the vicinity of the inner apparent horizon (again it is 
signaled by vanishing of the r v derivative), this horizon is displayed as a thick dash-dotted line. 
Note that the horizons "move" to approach asymptotically (v — * oo) each other. 




FIG. 3. The classical collapse. The radius r(u, v) as a function of the advanced time v along 
a sequence of u = const outgoing null rays. The retarded time u grows from the rightmost ray to 
the leftmost one. There is an outer apparent horizon located at about r ~ 16 (not shown) and the 
CH indicated by the vertical rays (see further discussion in Ref. |12,13[). 



To gain a further insight into the emerging spacetime we depict in Figure |] the radius 
contour lines in wf-plane. We show explicitly the inner apparent horizon (a thick dashed 
line) and the outer apparent horizon (a thick dot-dashed line). The outer horizon initially 
grows from zero, that corresponds to the formation of a black hole and then it shrinks 
signaling on the mass-evaporation. The inner horizon (CH) does not remain a contracting 
null hypersurface, as in a case without evaporation, but instead it expands outwards to meet 
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asymptotically the outer horizon. 

The spacetime of the EBH has a very different structure from the non-extremal one, in 
particular the outer and the inner horizons are located at the same radial coordinate. In 
the evolutionary context, where the initially non-extremal black hole evolves towards the 
EBH, the merge of the both horizons will occur in an infinite advanced time (as v — > oo). 
The numerical simulation can proceed only to finite values of v, and the horizons will only 
tend to approach each other. We indeed observe this during the simulation. The outer 
horizon (indicated by tick dotted line in Fig.|2|(b)) approaches the inner horizon (indicated 
by tick dashed line in Fig.Q(c)). The same picture is seen in Fig. |]. It should be noted 
that the shrinking of the outer apparent horizon is contained within a small variation of 
the u coordinates, as in a neutral collapse. That's because of a Kruskal-like nature of our 
coordinates. 
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FIG. 4. Left panel: lines of constant radius in uv-plane. The inner and the outer apparent 
horizons of an evaporating black hole are defined by vanishing of the r v derivative. The inner 
horizon (shown by thick dashed line) is approaching the outer horizon (the thick dash-dot line) in 
the asymptotically late advanced time, v — > do. Right panel: lines of constant radius in uv -plane 
in a vicinity of the event horizon. The outer apparent horizon is designated by the thick dash-dot 
line. This horizon contracts from about v ~ 50, indicating the mass-evaporation. On the both 
panels the radius is decreasing from the bottom to the top. 

It was interesting to compare the shrinking of the radius of the outer apparent horizon 
with the radius of this horizon, as it is calculated assuming the external RN metric. In the 
latter case the radius of the horizon is given by: 

r + = M+ \J M 2 - Q 2 . (9) 

Here and below M = Mar(u, v) and Q = Qah is the mass and the charge of the apparent 
black hole. In Figure § we depict the radius of the outer apparent horizon vs. the advanced 
time. One can observe that the numerically calculated radius of the horizon (designated 
by the circles in Figure) and the one obtained by the above formula (solid line) practically 
coincide. It must not be confusing, despite this remarkable result, the entire spacetime of 
the dynamically evaporating black hole is very different from the RN one. The radius of the 
inner horizon, denned in a RN spacetime as r_ = M — y/M 2 — Q 2 does not coincide with 
the calculated inner horizon. This is because of presence of energy influx which changes the 
mass-function inside the black hole. This mass-function is different from the external one, 
M AH . 
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FIG. 5. The location of the outer apparent horizon of the evaporating black hole as a function 
of the advanced time v. The horizon initially grows indicating the absorption of matter by the 
black hole. Subsequently (starting from v ~ 50) the horizon contracts, signaling on evaporation. 
The circles represent the numerically obtained radius of the horizon. The solid line represents the 
location of the horizon, which was calculated assuming RN metric, according to Eqn. (^). The 
notable similarity of the results shows that a charged black hole is indeed a RN black hole during 
its evolution, from a point of view of an observer at infinity. 

The Hawking-Bekenstein temperature of a RN black hole, as it measured in infinity is 
defined as: 



T 



VM> - Q 2 



(10) 



In Figure ^](a) we plot the dependence of the temperature of the dynamical black hole on 
the advanced time v. The temperature is seen to decrease. We note that the obtained black 
hole is close to a EBH as 1 - Q/M « 1/400. 

Another quantity of interest is the mass-loss rate. This can be calculated from a negative 
infalling energy flux (T vv ), which expected to drive black hole mass-evaporation according 
to: 

dM 



dv 



vvl \AH ■ 



(11) 



Figure |6|(b) depicts the explicitly calculated infalling energy flux (T vv ) along the outer ap- 
parent horizon as a function of an advanced time v. The negative influx approaches zero 
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value with the lapse of the advanced time v. This corresponds to the stable point in the 
evolution of the black hole. The black hole tends towards the EBH state, which would be 
reached at an infinite lapse of v. 




(a) ^ <b> V 

FIG. 6. (a)The Hawking-Bekenstein temperature T of a black hole as a function of an advanced 
time v along the outer apparent horizon. The temperature decreases as the black hole evaporates 
and approaches the EBH state. (b)The negative influx {T VV )\ AH as a function of v. The flux tends 
to the zero value in the late advanced time. 

B. Evaporation with Discharge 

In the next stage of the analysis we relax the assumption of a stable charge, allowing 
it to be radiated away via the Schwinger pair creation process. We assume that charged 
particles are not produced in the Hawking process. We examine analytically the fate of a 
simultaneously evaporating and discharging black hole. 

We suppose that the discharge begins only when the electric field of the collapsing matter 
approaches the critical value, E c , and there is no pair creation for the subcritical fields. This 
is justified since the rate of pair creation for subcritical fields is exponentially suppressed 
(see e.g. pf). 

Figure [I] is the QM-diagram. The black hole can only form if Q < M thus only the 
region below M = Q is allowed. M = Q is designated by the dashed line in this Figure. 
The solid line represents black holes for which: E = Q/r 2 ^ = E c along the outer horizon. 
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Below this line the electric field upon the outer horizon is subcritical and charged particles 



This line intersects the Q = M EBH line at Q = E~ l . After the intersection the line 
continues to the unphysical, Q > M, region. 

A trajectory of a black hole that evaporates, starting from an initial mass and charge, 
is represented by a horizontal line in this diagram. If this trajectory intersects the E = E c 
line it signals the onset of discharge. We depict such a trajectory (dash-dot line in the 
diagram) that all the three lines intersect at the same point, Q = E^ 1 , the point C in Figure 
[7|. The black holes that evolve along horizontal trajectories lying above the dash-dot line 
never experience a critical pair-creating field along their horizon. These black holes do not 
discharge but rather settle down to the stable EBH state. The fate of the black holes that 
do experience critical field is completely different. 

We show below that the following scenario holds. Black holes which evaporate along 
trajectories, lying below the dash-dot track reach the pair creating field at some moment. 
Then the black hole evolve along the E = E c line discharging and evaporating simultane- 
ously, decaying to a zero size. If the black hole evaporates exactly on the dash-dot track 
then it moves towards point C. Here we must recall that a small number of pairs will be 
produced even for a subcritical field. This small amount of pairs is sufficient to reduce (just 
a little) the charge of a black hole, pushing it away from the point C towards the decaying 
regime. This behavior will hold also for a small region above C as long as the pair creation 
is not suppressed exponentially. 



are not produced. Using @ we find the equation of this trajectory: 




(12) 
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1 2 

M 

FIG. 7. A QM-diagram. Black holes can exist only in region below Q = M track, designated 
here by the dotted line. The solid line in the diagram corresponds to the black holes which 
experience an extremal pair-creating electric field along their outer horizon. An evaporating black 
hole is presented in this diagram as a horizontal line. There is the limiting value of charge for an 
evaporating black hole, such that a black hole having charge above this value will not discharge 
but rather settle down to an EBH state. We depict a trajectory (dash-dot track) for a limiting 
black hole. 

Let us consider an evaporating black hole, moving along a horizontal line in the region 
below the dash-dot track. At some moment this black hole will approach the E(r + ) = E c = 
Q/r 2 c line. Suppose that the black have evolved slightly into the region above that line. The 
electric field along its outer horizon is now supercritical and pairs are intensively created. 
The created pairs reduce the charge of a black hole to the value Q c defined by Q c /r\ = E c . 
The pairs being created in the vicinity (1 — y- <C 1) of the outer horizon of the black 
hole are swallowed by it. This forces the increase, AM, of the mass of a black hole and 
therefore increase of the radius of the outer horizon r + — > r' + > r+. The electric field along 
r + is exactly critical and therefore along r' + is subcritical. As a result the black hole moves 
outwards away from the supercritical region. The black hole has evolved (M, Q) — > (M', Q c ) 
with M' = M + AM during the pair creating stage. In the next stage the black hole 
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evaporates until reaching (and maybe exceeding) again the critical field along its horizon. 
At the end of this phase the black hole would evolve (M',Q C ) — ► (M",Q C ). If the electric 
field along the outer horizon is supercritical the above story takes place once again and so 
on until the full decay The cycle (M, Q) — ► (M', Q c ) — ► (M", Q c ) is indicated by the cycle 
1-2-3 in the zoom-insert in Figure 0. The processes in nature are not discrete as we presented 
it but rather continuous and both phases take place simultaneously. 

The conclusion from the above considerations is clear - the black holes will evolve along 
the E(r + ) = E c track after approaching it. Along this track it satisfies 



dM 11 



+ 7 V^cQ . (13) 



dQ 4 V E C Q 4 
as obtained by differentiating Eq. ([12]). 

It is also of interest to calculate the mass- and the charge-loss rates of a black hole 

evolving along the E{r + ) = E c line. To address this question we implement the idea that if 

the electric field of a black hole exceeds the critical pair-creating value then the number of 

created pairs is not simply iVo = , but rather enormously amplified. This observation 

was previously made in |17j, tough the calculation we present below differs from that in 

Ref. |T7|] . The energy of the created pairs, e p , is extracted from the energy of the electric 



field of the black hole. If the electric field exceeds only slightly the critical value (that holds 
for 1 — y- <C 1 in our circumstance) then the pairs are created nearly at their rest energy 
with very little kinetic energy. This is in fact the definition of the critical pair-creating field. 
Therefore, the number of the created pairs can be estimated as N p ~ e p /m = AM/m and 
N p ^> No (holding for e/m ^> 1). The pairs are created in the region (r + ,r c ) with the 
proper thickness of 5r = r c — r + . The pairs constituting that region are in a state of the 
quasi-neutral plasma. Only the outermost and the innermost layers of this shell have an 



uncompensated charge of Q — Q c and the number of pairs in these layers is indeed A^ [17 . 

The rate of the pair creating process at a given radius within the pair-creating region 
can be approximated as 



2 



r =<" i-f^Ve-*. (14) 



From this we can estimate the time needed to create the N p pairs as 

iV 

Ar dis « ^ , (15) 

where V is the proper volume of the pair creating region. Using the above results we can 
estimate the mass increase and the charge loss by the black hole due to discharge. 

The charge loss is estimated as Q = ~ ~ ^J^ C TV. The total energy of pairs 
and therefore the mass increase, AM, of the black hole is estimated knowing the initial, 
Q/r 2 , and the final, Q c /r 2 electric field configurations 

O 2 - O 2 

AM « — ^ . (16) 

2r + 

Here we have ignored the energy carried by the N particles, repulsed to infinity since 
N p ^> N . Using Eq. (|16"D to calculate N p = AM/m one obtains, using 5r <C r + , the charge 
loss 

<>•* < 17 > 

In the evaporation phase the mass is radiated away by the Hawking process. In order 
to estimate this rate of mass loss we assume that the black hole radiates as a black body 
with temperature T, defined by QTUD, and surface area A = Airr 2 ^. Hence we obtain the 
approximate mass-evaporation rate: 

M ev = H w aT 4 A , (18) 

where a is the Stefan-Boltzmann constant. 

Both processes of evaporation and discharge are simultaneous. The evaporating black 
hole approaches and exceeds slightly, say by Sr, the critical radius r c . The value of Sr, 
which defines the volume of the pair-creating region according to V — Anr^dr, will be 
determined by the requirement that the ratio M/Q is given by (0). That is the black 
hole evolves nearly along the E(r + ) = E c line. The rate of the mass loss is given by 
M = M ev - M dis = H - w H - mTV. The rate of the charge loss is given by Eq. flT7|). 
Therefore, from 
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\dQ ) Q (^)mTV 
assuming that (4tt) = (4?f) one obtains 

to V dQ ) \dQ J E{r+)=E C 



dM\ M H-mTV , s 

(19) 



M = H 



• i 



(20) 



r+(dM_\ | i 

\ Q \ d Q J E(r + )=E C J 

With the use of Eqs.(|l"2"D, fll3|) and fll8|) this determines the rate of the mass loss in the 
combined discharge-evaporation process. We have also to perform the necessary consistency 
check which assures that Sr/r + <C 1. From Eqs. (|T9|) , fl20|) and ( |I4D with T(E) « T(E C ) one 
finds that Sr/r + 1 if 1/Mm <C 1 i.e. if X/r + 1, where X is the Compton wavelength of 
the created particles. The last inequality holds until the very late stages of the decay when 
the discharge occurs also via the Hawking process. Hence until the late phases of decay our 
calculation remains reliable. 

Since £ (%) is positive and both £ > 1 and (%\ > 1, the rate (1201) 

Q WO/E(r + )=B c <2 \ d Q J E{r+)=E C ' 

of decay of charged black hole is slower than that of the uncharged evaporating one M < H 
and together H > M > H/2. Therefore the life-time of the evaporating and discharging 
black hole is longer than the life-time of an evaporating neutral black hole of the same mass. 
This is readily understood - what we have after all is the conversion of electrostatic energy 
of a black hole to the mass, which is radiated away as a Hawking radiation. Since that 
effective mass exceeds M (approximately by the energy of the electric field), the process of 
decay lasts more time. 

Clearly the decay as we describe it is valid only above the Planck scales, where our 
(semiclassical) approximation breaks down in favor of full Quantum gravity. 



V. SUMMARY AND CONCLUSION 

We have investigated the formation and evaporation of charged black holes. In the first 
stage of analysis we assumed stability of the charge, while allowing an evaporation and a 
decrease of the mass via Hawking radiation. The evaporation was modeled by introducing 
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the expectation value of a 2D stress-energy tensor of a 2D quantized scalar field as a source to 
the semiclassical Einstein equations. By doing so, we disregard possible contributions which 
arise from the dimensional reduction procedure ||. The source which we utilized is (at least) 
a part of the full 4D stress-energy tensor, and it gives a feeling of the full back-reaction. 

We find that the evaporating black hole loses its mass and approaches asymptotically a 
stable endpoint corresponding to the extremal black hole. Starting from an initially regular 
matter distribution we have dynamically obtained a non extremal black hole and than, after 
evaporation, a near extremal one. The outcome of our analysis is the entire^ dynamical 
spacetime of such a black hole. 

The spacetime structure of a non extremal and an extremal black holes are very different. 
The most striking feature of the latter is that its inner and outer apparent horizons are 
located at the same radial coordinate. In the dynamical context we observe that the horizons 
tend to approach each other in a late advanced time. The actual merge will occur indeed 
in an infinite (advanced) time, i.e. as v —>■ oo. Thisis in agreement with the unattainability 
formulation of the third law of the black hole thermodynamics. The stability of the final 
EBH state is manifest from the vanishing of the negative energy influx, (T vv ), governing 
the black hole evaporation according to Eq. QTTJ) . The Hawking-Bekenstein temperature 
decreases as the black hole approaches the EBH state. 

We have also considered a collapse with evaporation accompanied by a discharge due 
to instability of the vacuum in strong electric fields. In this case there is a limiting value 
of charge for a given mass such that black holes with a charge above this value would not 
discharge but rather evaporate towards a stable EBH state. Black holes with a charge 



2 Actually, our solution does not include the origin of coordinates, r = 0, and the asymptotically 
late advanced time region. The first is not included in a sake of simplicity of the numerics, while the 
second is not included since it has an infinite extension. In any case the region which is obtained 
supply enough information to draw the conclusions, which we have drawn. 
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below the limiting charge evaporate and discharge simultaneously, decaying to zero size. 
It is interesting to note that such black holes are pushed away from the EBH state. The 
discharge thus prevents from small black holes to become extremal. The life time of an 
evaporating and discharging black hole is longer then that of an uncharged black hole of the 
same mass that decays to zero emitting Hawking radiation. 

Note added: When this work was completed its brought to our attention that there are 
related works (T^j considering evaporation of 2D RN black holes. 
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